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В статье рассматривается подход к решению проблемы снижения пропусков в обнаружении изменений 
наземных объектов по данным спутникового дистанционного зондирования Земли с автоматизацией 
соответствующего процесса в системах детерминированного тематического мониторинга. 
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Введение 


Одной из тенденций в создании систем детерминированного тематического 
мониторинга, использующих данные спутникового дистанционного зондирования 
Земли (ДЗЗ), является обеспечение повторяющегося от космической съемки к съемке 
объема обработки данных. Подчеркивая отличие такого рода систем от других, их 
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называют промышленными [1], основываясь на использовании в их построении пол- 
ностью автоматизируемых (программно реализуемых) методов. В других случаях их 
именуют системами постоянного и даже — оперативного мониторинга. 

Можно было бы класс таких систем назвать и ведомственным, подразумевая, 
что системы постоянного контроля определенных участков поверхности Земли должны 
эксплуатироваться в интересах министерств сельского хозяйства, архитектуры и стро- 
ительства, чрезвычайных ситуаций, обороны и др. 

Выделяя при этом значимость тенденции полной автоматизации для систем 
мониторинга любого класса, можно увидеть, что для интересующих систем она ока- 
зывается неотвратимой. В силу повторяемости объема и последовательности решаемых 
задач последние должны «обеспечить качественный анализ данных без участия чело- 
века» [1]. Неотвратимость автоматизации вызывается рутинным характером обработки 
регулярно получаемых значительных объемов спутниковой информации разновре- 
менных съемок для больших территорий. Оставлять при этом операторское участие — 
значит вносить в результаты тематического мониторинга субъективность, снижающую 
уверенность в получении объективных и сопоставимых разновременных данных, 
соглашаться с недопустимо высокой стоимостью эксплуатации. 

Все это говорит в пользу такой наиболее удачной общей характеристики и наз- 
вания рассматриваемых систем как «системы детерминированного тематического 
мониторинга (СДТМ)». По соответствующему характеру решаемых задач им должно 
сопутствовать достижению главной цели — автоматическое обнаружение возникающих 
изменений, что предполагает построение СДТМ на основе использования данных 
разновременных космических съемок контролируемого участка Земли. 

Достижение рассмотренных качеств СДТМ должно начинаться с автоматизации 
выделения и обработки заданной области мониторинга из состава полосы захвата 
отдельной спутниковой съемки земной поверхности, на что и ориентирована настоящая 
статья. 


Область мониторинга в космических съемках 
земной поверхности 


Обеспечение автоматического характера СДТМ требует нового подхода к раз- 
мерам изображений объектов съемки, подвергаемых обработке в их составе. Причиной 
этого обычно является зависимость потребностей целого ряда тематических систем в 
охвате съемкой больших участков земной поверхности при высокой разрешающей 
способности. В результате этого интересующая площадь может быть охвачена съемкой 
только в ряде пролетов спутника. 

Точно так же, если даже область мониторинга не превышает размеров изобра- 
жения одной съемки, она может быть получена в последовательных пролетах. Причем 
на очередном витке спутника в силу случайных факторов не гарантируется, что участок 
наземных интересов полностью войдет в зону захвата. Наряду с этим при полном 
или частичном вхождении соответствующее ему изображение выполненной съемки 
будет содержать второстепенные, не интересные для мониторинга участки. 

При получении изображения области наземных интересов по частям, то есть по 
разным съемкам, приходится принимать во внимание случайность размеров этих 
частей. Она оказывается одной из причин использования в полной сборке области 
мониторинга квалифицированного оператора, что противоречит рассмотренным прин- 
ципам построения СДТМ. 
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Так, часть изображения ©, интересующей области 62, полученная в эталонной 
съемке, не может играть роли целой (<. = ©). В то же время часть этой же области <>,, 
но полученная в очередной съемке с некоторым перекрытием ее (ДО. ) с областью ©, 
в качестве дополнения, нуждается в удалении из ее состава этого перекрытия: 


О) =О, /АО.,. Размеры АО», случайны от сеанса к сеансу съемки. Всю интере- 
сующую область можно охватить в результате интерактивных действий с рассмот- 
ренной парой съемок, выполняя объединение: @ =О 0). 

Рассмотренные причины и определяют в качестве первоочередной задачи авто- 
матизацию формирования области наземных интересов, называемой областью мони- 
торинга. В то же время главная решаемая задача системы — обнаружение изменений 


состояния этой области — не позволяет обойти оценку его возможного влияния на 
представления о ней. 


Влияние решения основной задачи СДТМ 
на размеры области мониторинга 


В соответствии со своим назначением главной задачей СДТМ является обнару- 
жение возникающих изменений, что естественно предполагает использование данных 
разновременных космических съемок контролируемого участка Земли. Одна из них 
доставляет опорное (эталонное) изображение, а вторая — изображение текущей съемки. 
Решение задачи обнаружения представляет собой поиск изменений состояния отоб- 
ражаемых объектов текущей съемки путем сравнения очередного изображения с 
предшествующим ему по времени съемки (эталонным). 

Наиболее простым и позволяющим надеяться на автоматизацию обнаружения 
является метод разностных изображений. Для его реализации наиболее целесообразно 
воспользоваться представлением матрицы разности изображений разновременной 
съемки в виде взаимодополняющей пары изображений: 


АКх, у), если АЦх, у) >0; 


ДЕ, (х, у) = 
0 востальных случаях; 


АЦх, у), если АЦх, у) <0; 


0 востальных — случаях, 


ДЕ, (х,у)= (1) 
где АМх, у) = Цх, у) - [.(х, у) — разностное изображение; 
[х, у), [-(х,Уу) — текущее и эталонное изображения как матрицы яркостей 
пикселей, имеющих пространственную привязку (х, у) в их составе; 
АЕ, (х,у) — изображение, отражающее возрастание яркостей при изменениях; 


АЁ,(х, у) — изображение, отражающее уменьшение яркостей при изменениях. 


На рис. 1 представлены изображения разновременной съемки одного и того же 
участка земной поверхности в районе г. Брест размером 9,8х9,3 км, выполненной 


спутником ОшскВИА (май, июль) и соответствующая им взаимодополняющая пара, 
представляющая разностное изображение. 
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в) 


Рисунок 1 — Изображения разновременной космической съемки 
(а май, б — июль) и взаимодополняющая пара разностных изображений 
(в — увеличения яркостей, г — уменьшения яркостей) 


С учетом введенной выше терминологии приведенное изображение более ранней 
(майской) съемки (рис. 1а) принципиально могло бы играть роль эталона для входящей 
в его состав области мониторинга или ее части. Изображение съемки, выполненной в 
июле (рис. 16), следовало бы отнести при этом к текущей съемке того же участка. 

Рис. | позволяет обратить внимание: 

— на реально имеющую множественность изменений, наблюдаемых в составе 
одного и другого взаимодополняющих изображений; 

— на одновременное наличие крупноразмерных и мелких яркостных изменений; 

— на близость к фоновым яркостям отдельных изменений; 

— на яркостные отличия между собой компонент изменений в пределах изобра- 
жений каждого знака яркостей. 

С точки зрения выполнения задач мониторинга, предполагающих реакцию на 
изменения, каждое изображение рассмотренной пары может представлять самостоя- 
тельный интерес. При этом обнаружение любого изменения обычно влечет за собой 
необходимость его раздельной сегментации, определения площади, диапазона яркостей 
и координат. 

Одновременно с обеспечением зрительного наблюдения особенностей изменений 
в пределах изображения космической съемки рис. 1 позволяет продемонстрировать 
упомянутую в п. 1 возможность несовпадения области мониторинга со всем изобра- 
жением космической съемки. В частности согласно рис. Та, б в качестве такой области 
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может быть принят, например, только бассейн реки Лесная (наиболее хорошо зри- 
тельно выделяется на рис. 16), характеризуемой на снимке направлением течения с 
юго-запада на северо-восток с поворотом почти на север в центре отображаемого 
района. Точно также это мог быть и один из лесных массивов или их обособленные 
участки. Причем в любом из таких случаев изменения, обнаруживаемые вне выбранных 
областей изображения съемки, оказывались бы не интересующими наблюдателя. 

Отсюда для контроля самих областей, казалось бы, необходимо использование 
масок, выделяющих эти области из состава района съемки. Но сложный характер 
границ природных образований приводил бы к геометрической сложности и к сопут- 
ствующим трудностям использования масок для непосредственного сопоставления 
их разновременных изображений при выделении изменений. Кроме того, решение 
такой задачи может сопровождаться упомянутыми (п. 1) проблемами последовательной 
съемки области мониторинга по частям. 

Однако для полноты понимания взаимосвязи выбора области мониторинга и 
обеспечения автоматического обнаружения изменений в ее пределах необходимо 
рассмотреть возможность решения соответствующей задачи применительно к условиям, 
представленным изображениями рис. 1. 

Наиболее очевидным подходом к обнаружению яркостных изменений является 
комплексирование разностного изображения области мониторинга с пороговой обра- 
боткой [2], [3]. Методы пороговой обработки в силу простоты и соответствия зрительно 
достигаемым результатам наиболее очевидны. 

Например, в случае одной компоненты изменений (©„) знание яркостного 


порога АЁ,,„ для любого изображения взаимодополняющей пары обеспечивает 


выделение соответствующих ему яркостей: 
АЕ, «(х, у), если АГ, „(х, У) > АР, д; 


АЁ, а (х, У ЕО,„) = 
0, в остальных случаях. 


(2) 

Определение самого порога предполагает построение гистограмм по каждому 
изображению изменений рассматриваемой взаимодополняющей пары, как это для 
изображений рис. | приведено на рисунке 2. 

Каждая гистограмма при этом отражает как совпадающие части изображений 
разновременной съемки, так и части, подвергнутые изменениям. При этом фон обра- 
зуется как в результате полного совпадения неизмененных участков, так и воздействия 
яркостного шума съемок. Содержащий нулевые и близкие к нему яркости в составе 
гистограммы он формирует группировку пикселов, порожденных распределениями 
яркостного шума и других источников, у оси ординат гистограммы. 


6000000 5000000 
5000000 4000000 
4790080 3000000 
3000000 

0066 2000000 
1000000 1000000 
0 т 0 

ме о [>] со г (К) [Ге] < ее (А) со г >] те) 

сб Го) со + ыи г [52] [Ге] [> и, ыя гг 

< < < < + - 

а) 6) 


Рисунок 2-— Гистограммы взаимодополняющей пары разностных изображений 
(рисунка 1 В, Г): а - гистограмма составляющей изменений, вызванных возрастанием 
яркостей; б — гистог рамма изменений, вызванных их уменьшением 
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В отличие от этого яркостные изменения в составе гистограммы представляются 
в виде выборок пикселов, образованных некоторыми распределениями вероятностей 
отражающих свойств подстилающих поверхностей изменений в диапазонах их яркостей, 
соответствующих выражению (2). Отсюда использование гистограмм разностного 
изображения в определении порога должно сводиться к оценке яркости, разделяющей 
в составе гистограммы фон и измененную часть. Зрительно это выполняется доста- 
точно просто, а для автоматического решения известно большое число различных 
методов [4]. 

Существующее многообразие объектов съемки по количеству компонент и по 
их физическим особенностям можно считать случайными факторами. Их влияние 
определяет случайный характер гистограмм взаимодополняющих пар разностных 
изображений. 

Главная особенность гистограмм, обнаруживаемая на рис. 2, относится к харак- 
теру составляющей рис. 2 а. Здесь — это отсутствие явно наблюдаемых изменений в 
составе упомянутой гистограммы, зрительно отмечаемых на рис. 1 в. Соответствую- 
щими причинами такого явления зрительный анализ исходного изображения (рис. | в) 
заставляет считать: 

— слабое отличие изменений от фона; 

— равномерное яркостное заполнение имеющимся множеством яркостных изме- 
нений всего интервала их, примыкающего к фону. 

Сопутствующими последствиями указанного отсутствия является пропуск име- 
ющих место изменений. 


В отличие от этого гистограмма составляющей уменьшения яркостей Н(АЁ,) 


(рис. 2 6) демонстрирует одну компоненту изменения, отделенную от фоновой части 
небольшой седловиной. Соответствующий этому выделению изменений порог сос- 


тавляет АЁ и = 73. В то же время эта компонента, судя по зрительным оценкам рис. 1 г, 


объединяет в себе все множество наиболее ярких компонент. претерпевших изменения. 
Менее яркие изменения могут так же, как и в случае рис. 2 а, поглощаться фоном. 
Для подтверждения этого достаточно выполнить сегментацию, основанную на бина- 
ризации изменений, превышающих порог: 


1 ири Ах, у) > АЁк; 
В(х, у) = 0 (3) 
в остальных случаях. 

Соответственно приведенному выше значению порога на рис. 3 представлена 
сегментация изменений, связанных с уменьшением яркостей в изображении очередной 
съемки по сравнению с предшествующим ему по времени изображением эталона. 
Причина того, почему для изображения возрастания яркостей такая бинаризация 
невозможна, состоит в отсутствии явно выраженных изменений в составе соответст- 
вующей гистограммы (рис. 2 а). Согласно рис. 3 очевидно, что большая площадь 
области мониторинга (здесь она предполагалась равной всей площади, охваченной 
снимком) в гистограммном представлении разностного изображения, является при- 
чиной потери индивидуальности отдельных изменений. 

Убедиться в достижении большей детальности для тех же изображений по гис- 
тограммам взаимодополняющей пары (рис. 1 в, ги 1 а, 6) достаточно просто, рассмо- 
трев особенности малоразмерных фрагментов. 
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Рисунок 3 — Сегментация изменений (6) разностного изображения (а), 
отображающего уменьшение яркостей изображений очередной съемки 
по сравнению с изображением эталона 


Так, на рис. 4 а, 6 приведен фрагмент разновременной съемки размером 2х1 км 
из состава изображений рис. 2 а, б (северо-восток территории). Рис. 4 в, г отображают 
при этом взаимодополняющую пару разностного изображения фрагмента, рис. 4 д, е — 


соответствующие им гистограммы Н(АГ,), Н(АГ,), а рис. 4 ж,зи 4 и - результаты 


сегментации изменений возрастания и уменьшения. 

Рассмотрение приведенных данных по фрагменту (рис. 4) в отличие от гисто- 
грамм разностного изображения всей области съемки (рис. 2) и сегментации ее изме- 
нений свидетельствует о достаточно четком отображении детальных изменений, как 
с возрастанием, так и с убыванием яркостей. Таким образом, с точки зрения избав- 
ления от пропусков в обнаружении отдельных изменений, характерных для большой 
области мониторинга, оправдывается необходимость уменьшения ее размеров. 


Дискретизация области мониторинга 


Несмотря на достигнутый эффект, обнаруженная тенденция уменьшения области 
мониторинга вступает в противоречие с практически возникающей потребностью 
максимального использования возможностей космической съемки по охвату ею земной 
поверхности. 

Выходом из рассмотренной ситуации, обеспечивающей снижение пропусков в 
обнаружении изменений и одновременное решение задачи автоматического форми- 
рования и обработки области мониторинга, является ее дискретное представление. 

В качестве фрагментов дискретизации легко предложить использование мно- 
жества изображений равного размера, одинаково ориентированных по сторонам света, 
но приобретающих индивидуальную топографическую привязку в области монито- 
ринга. При этом, чем меньше площадь фрагмента, тем со всей очевидностью более 
точным может быть описание такой области контроля изменений. Поэтому любой 
участок наземных интересов всегда может быть построен из расположенных вплотную 
малых по размерам фрагментов анализа, как это показано на рис. 5. 
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Рисунок 4 — Гистограммы фрагмента 2х1 км взаимодополняющей пары 
разностных изображений: эталон, текущее (а, 6); изменения возрастания, изменения 
уменьшение яркостей (в, г) и соответствующие им гистограммы (д, е); сегменты 
изменений с возрастанием яркостей (ж, з) и сегмент изменений — с уменьшением (и) 
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Рисунок 5 — Примеры дискретизации области мониторинга на реальном 
изображении спутниковой съемки территорию Беларуси 


Убедительность вывода о повышении детальности контроля изменений, снижение 
числа пропусков с уменьшением площади дискрета-фрагмента подтверждается экс- 
периментально при обращении к конкретным методам обнаружения с учетом каждого 
тематического содержания мониторинга (лес, водная поверхность, с/х посевы и т.п.). 
В то же время после выбора фрагмента дискретное формирование всей области мони- 
торинга для конкретной СДТМ не требует коррекции от случая к случаю. Это одно- 
кратно выполняемый для нее элемент настройки. 

Основой для такой настройки должна быть замена карты местности, используемой 
в СДТМ при установке области мониторинга, ее фотопланом [5]. Применение фото- 
плана, привязанного к местности и представляющего ее плановое панхроматическое 
изображение с той же разрешающей способностью, что и текущие космические съемки, 
может быть легко поддержано программными средствами, облегчающими расстановку 
фрагментов (выбор положений и осуществление перемещений). В результате каждый 
фрагмент должен приобретать свою топографическую привязку (координаты). При этом 
ввод очередного участка фотоплана (съемка по частям) достаточно просто сопро- 
вождается отображением всех расставленных до этого фрагментов, попавших в зону 
соответствующей съемки. В результате появляется возможность расстановки новых 
фрагментов на неохваченной до этого части области мониторинга. 

С учетом выполненной дискретизации и описания области мониторинга подго- 
товка к обнаружению изменений должна преследовать две цели: 

— подготовительная, состоящая в выделении и импортировании на хранение в 
качестве эталонных (опорных) изображений всех тех фрагментов, которые составляют 
область мониторинга и отражают исходное или целесообразное состояние отобража- 
емых участков земной поверхности; 

— функциональная, состоящая в выделении и импортировании на обнаружение 
изменений текущих фрагментов области мониторинга, для которых подтверждается 
наличие эталонных изображений. 

Кроме ориентации на детальный контроль изменений и рассмотренную простоту 
поэлементного обнаружения всей области или части ее в пределах очередной съемки 
района, использование ее дискретизации в силу равенства размеров отдельных фраг- 
ментов приводит к унификации рекуррентно повторяющейся обработки. 
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Выводы 


1. Выяснение требуемых качеств системы детерминированного тематического 
мониторинга (СДТМ), состоящих в повторяемости объема и последовательности 
решаемых задач обработки, исключении участия операторов из этого процесса, пока- 
зало, что автоматизация сопутствующих процессов должна начинаться с автоматизации 
выделения и обработки заданной области мониторинга из состава полосы захвата 
отдельной спутниковой съемки земной поверхности. 

2. Рассмотрение особенностей обнаружения возникающих изменений, как главной 
задачи СДТМ, позволило установить, что большие по размерам области мониторинга 
могут характеризоваться отсутствием явно наблюдаемых множественных изменений 
в составе гистограмм, участвующих в реализации пороговой обработки изображений. 
Причинами этого являются пропуски в обнаружении компонент изменений, диапазон 
яркостей которых не отличается или отличается незначительно от яркостного диапа- 
зона фона, образованного помехами регистрации и ошибками геометрической и яркос- 
тной коррекции изображений разновременной съемки. Снижение числа пропусков 
обеспечивает тенденция уменьшения размеров области мониторинга. 

3. Стремление использования преимуществ в больших охватах земной поверх- 
ности спутниковой съемкой вступает в противоречие с установленной тенденцией 
уменьшения размеров области мониторинга. Устранение противоречия с одновре- 
менным повышением качества аппроксимации области мониторинга и снижением 
вероятности пропусков в обнаружении изменений достигается предлагаемой дискре- 
тизацией этой области, принципов выбора фрагмента дискретизации и особенностей 
их системной реализации в достижении целей. 
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